Binding-induced functional-domain motions in the Argonaute characterized by adaptive advanced sampling.

Argonaute proteins in combination with short microRNA (miRNAs) can target mRNA molecules for translation inhibition or degradation and play a key role in many regulatory processes. The miRNAs act as guide RNAs that associate with Argonaute and the complementary mRNA target region. The complex formation results in activation of Argonaute and specific cleavage of the target mRNA. Both the binding and activation processes involve essential domain rearrangements of functional importance. For the Thermus Thermophilus Argonaute (TtAgo) system guide-bound (binary) and guide/target-bound (ternary) complexes are known but how the binding of guide and target mediate domain movements is still not understood. We have studied the Argonaute domain motion in apo and guide/target bound states using Molecular Dynamics simulations and a Hamiltonian replica exchange (H-REMD) method that employs a specific biasing potential to accelerate domain motions. The H-REMD technique indicates sampling of a much broader distribution of domain arrangements both in the apo as well as binary and ternary complexes compared to regular MD simulations. In the apo state domain arrangements corresponding to more compact (closed) states are mainly sampled which undergo an opening upon guide and guide/target binding. Whereas only limited overlap in domain geometry between apo and bound states was found, a larger similarity in the domain distribution is observed for the simulations of binary and ternary complexes. Comparative simulations on ternary complexes with 15 or 16 base pairs (bp) formed between guide and target strands (instead of 14) resulted in dissociation of the 3'-guide strand from the PAZ domain and domain rearrangement. This agrees with the experimental observation that guide-target pairing beyond 14 bps is required for activation and gives a mechanistic explanation for the experimentally observed activation process.

Influence of a Ser111-phosphorylation on Rab1b GTPase conformational dynamics studied by advanced sampling simulations.

Rab GTPases constitute the largest branch of the Ras protein superfamily that regulate intra-cellular membrane trafficking. Their signaling activity is mediated by the transition between an active GTP-bound state and an inactive GDP-bound state. In the inactive state the switch I and II segments adopt largely disordered flexible conformations, whereas in the active state these regions are in well-defined conformations. The switch I and II states are central for recognition of Rab GTPases by interacting partners. Phosphorylation of the Rab1b-GTPase at residue Ser111 (pS111) results in modulation of the signaling activity due to alterations of the protein interaction interface and also due to modulation of the conformational flexibility. We have studied the flexibility of native and pS111-Rab1b in complex with GTP or GDP using extensive Molecular Dynamics (MD) simulations and an advanced sampling method called DIhedral Angle-biasing potential Replica-Exchange Molecular dynamics (DIA-REMD). The DIA-REMD method promotes backbone and side chain dihedral transitions along a series of replica simulations in selected protein segments and through exchanges also improves sampling in an unbiased reference simulation. Application to the Rab1b system results in significantly enhanced sampling of different switch I/II conformational states in the GDP-bound Rab1b state. The pS111 modification is found to reduce the conformational flexibility even in the presence of GDP, which may influence signaling activities. The stabilizing effect can be attributed to the formation of additional surface salt bridges between Arg-residues and pS111 not present in the native structure. The DIA-REMD method could be a valuable approach for studying also other signaling proteins that contain flexible segments.

Specificity of AMPylation of the human chaperone BiP is mediated by TPR motifs of FICD.

To adapt to fluctuating protein folding loads in the endoplasmic reticulum (ER), the Hsp70 chaperone BiP is reversibly modified with adenosine monophosphate (AMP) by the ER-resident Fic-enzyme FICD/HYPE. The structural basis for BiP binding and AMPylation by FICD has remained elusive due to the transient nature of the enzyme-substrate-complex. Here, we use thiol-reactive derivatives of the cosubstrate adenosine triphosphate (ATP) to covalently stabilize the transient FICD:BiP complex and determine its crystal structure. The complex reveals that the TPR-motifs of FICD bind specifically to the conserved hydrophobic linker of BiP and thus mediate specificity for the domain-docked conformation of BiP. Furthermore, we show that both AMPylation and deAMPylation of BiP are not directly regulated by the presence of unfolded proteins. Together, combining chemical biology, crystallography and biochemistry, our study provides structural insights into a key regulatory mechanism that safeguards ER homeostasis.

Proteolysis of Rab32 by Salmonella GtgE induces an inactive GTPase conformation.

Rab GTPases are central regulators of intracellular vesicular trafficking. They are frequently targeted by bacterial pathogens through post-translational modifications. Salmonella typhimurium secretes the cysteine protease GtgE during infection, leading to a regioselective proteolytic cleavage of the regulatory switch I loop in the small GTPases of the Rab32 subfamily. Here, using a combination of biochemical methods, molecular dynamics simulations, NMR spectroscopy, and single-pair Förster resonance energy transfer, we demonstrate that the cleavage of Rab32 causes a local increase of conformational flexibility in both switch regions. Cleaved Rab32 maintains its ability to interact with the GDP dissociation inhibitor (GDI). Interestingly, the Rab32 cleavage enables GDI binding also with an active GTP-bound Rab32 in vitro. Furthermore, the Rab32 proteolysis provokes disturbance in the interaction with its downstream effector VARP. Thus, the proteolysis of Rab32 is not a globally degradative mechanism but affects various biochemical and structural properties of the GTPase in a diverse manner.

Phosphorylation of Ser111 in Rab8a Modulates Rabin8-Dependent Activation by Perturbation of Side Chain Interaction Networks.

GTPases are key players during cellular signaling. Phosphorylation of Rab proteins, which belong to the Ras superfamily of small GTPases regulating intracellular transport, has been implicated in the pathogenesis of Parkinson's disease. For Rab8a, it was shown that serine 111 phosphorylation (pS111) is dependent on the protein kinase PINK1 and that mimicking the phosphorylation at S111 by a serine/glutamate substitution (S111E) impaired Rab8a activation by its cognate nucleotide exchange factor (GEF) Rabin8. However, Ser111 is not part of the interface of the Rab8a:Rabin8 complex. Here, we performed comparative molecular dynamics and free energy simulations on Rab8a and Rab8a:Rabin8 complexes to elucidate the molecular details of how pS111 and S111E may influence the interaction with Rabin8. The simulations indicate that S111E and pS111 establish an intramolecular interaction with arginine 79 (R79). The interaction persists in the complex and perturbs a favorable intermolecular salt-bridge contact between R79 in Rab8a and aspartate 187 in Rabin8. Binding free energy analysis reveals that S111E and pS111, as well as the R79A mutation, drastically decrease the binding affinity for Rabin8. Combining the R79A mutation with S111E or pS111 nearly diminishes Rab8a-Rabin8 binding. In vitro experiments confirm our computational results showing a >80% decrease in the nucleotide exchange rate of the respective Rab8a mutants in the presence of Rabin8 compared to that of the wild type. In addition to insights into how S111 phosphorylation of Rab8a influences GEF-mediated activation, the simulations demonstrate how side chain modifications in general can allosterically influence the surface side chain interaction network between binding partners.